Ferroelectrically induced weak-ferromagnetism in a single-phase multiferroic by design 
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We present a strategy to design structures for which a polar lattice distortion induces weak ferro- 
magnetism. We identify a large class of multiferroic oxides as potential realizations and use density- 
functional theory to screen several promising candidates. By elucidating the interplay between the 
polarization and the Dzyaloshinskii-Moriya vector, we show how the direction of the magnetization 
can be switched between 180° symmetry equivalent states with an applied electric field. 
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The rational design of new materials with emergent 
properties is a riveting challenge today in materials 
physics. It begins with understanding a mechanism to 
control the interplay between diverse microscopic degrees 
of freedom in order to create targeted macroscopic phe- 
nomena and ends with the discovery or design of new ma- 
terial realizations. When combined with first-principles 
density-functional theory, this approach provides an ef- 
ficient strategy to survey the vast space of possible ma- 
terials to target for synthesis. For example, new multi- 
ferroics in which magnetism coexists with ferroelectricity 
have been discovered where magnetic order itself induces 
ferroelectricity. Through this specific spin-lattice interac- 
tion, it is readily possible to control the direction of the 
electrical polarization with a magnetic-field [H, l2|. An 
equally fundamental but technologically more relevant 
problem that has received far less study is the electric- 
field control of magnetism 0, 0, H, 0]. In particular, 
the electric-field switching of a magnetization between 
180° symmetry equivalent states has yet to be demon- 
strated. The most promising direction to achieve this in 
single-phase materials involves a ferroelectric distortion 
inducing weak-ferromagnetism [7',[i,i,[lQ]. Discovering 
a prototypical structure for which this approach might 
be realized, however, has remained elusive. 

Weak-ferromagnetism (wFM) is the phenomenon 
whereby the predominantly collinear spins of an an- 
tiferromagnet cant in such a way as to produce a 
residual magnetization (M). It can arise as a rela- 
tivistic correction to Anderson's superexchange, i.e., 
the Dzyaloshinskii-Moriya (DM) interaction [ill, [lH, 
£^DM=^ij-SiXSj where D is the Dzyaloshinskii vector. 
A ferroelectric (FE) distortion can induce wFM when 
the phenomenological invariant £^plm'^P-(LxM) is al- 
lowed in the energy of the antiferromagnetic-paraelectric 
(AFM-PE) phase, where P and L are the polarization 
and AFM vector respectively. Due to this term a cou- 
pling between the sign of P and M, for fixed L, is evident. 
It is important to realize that a FE can still exhibit wFM 
without the FE distortion per se causing the wFM, i.e., 
without E'pLM in the corresponding PE phase. For exam- 
ple, although BiFeOa (the most widely studied multifer- 
roic) can display wFM, such an invariant does not exist 



as previously shown from explicit first-principles calcu- 
lations [13] and as we argue below from symmetry. The 
challenge then is understanding how to start with the 
microscopic properties of £^dm and build a material cap- 
turing the macroscopic physics of £^plm- In this Letter 
we present for the first time design criteria that facilitate 
this process. In general the criteria target a structure 
for which wFM is symmetry-forbidden in the PE phase 
but symmetry- allowed in the FE phase [4, J^]. Below we 
use this to identify a class of multiferroic oxides as po- 
tential realizations and use density-functional theory to 
screen several promising candidates. One complication 
arises in magnetic systems with broken spatial inversion 
such as a FE: a symmetry allowed interaction [l3| niay 
in some cases generate a long- wavelength magnetic spiral 
canceling the net M, e.g., as in (bulk) BiFeOa. Here we 
consider the situation where this inhomogeneous state is 
suppressed, e.g., when single- ion anisotropy is large [l5[. 

Our strategy is to formulate the problem in terms of 
a structural-chemical criterion and a magnetic criterion. 
To illustrate the idea we consider a two-sublattice antifer- 
romagnet such as BiFeOs. A synopsis of the structural- 
chemical criterion is as follows: start with a PE structure, 
decorate the lattice with spins such that the midpoint be- 
tween two spins coincides with a site having inversion (X) 
symmetry (so that D = {)hy symmetry, i.e., Moriya's first 
rule [12]), place a FE-active ion at anX-site. Compounds 
which satisfy this criterion are quite intriguing because if 
there are no other symmetry elements that would forbid 
wFM - the magnetic criterion - all that is required to 
induce a non-zero D and wFM is to remove the X center 
by controlling the off-centering of the FE-active ion either 
by temperature, pressure, or electric-field. The magnetic 
criterion is primarily a question of how the spins order 
and of their direction with respect to the crystallographic 
axes. Notice our microscopic based criteria implies that 
L is odd under X, which is precisely what is required con- 
sidering macroscopic phenomenology, i.e., £^plm- These 
criteria facilitate evaluating known compounds or when 
combined with crystal chemistry principles designing new 
prototypes. For example, using PE BiFeOs as a starting 
point, we apply them to design a novel structure. 

In PE BiFeOs, crystallographic space group R3c, the 
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FIG. 1: Paraelectric ABO3 compounds, space group, R3c, 
BiFeOs and FeTiOs. A-site: WyckofF position 2a, site sym- 
metry 32. B-site: Wyckoff position 2b site symmetry, 3. 



Bi-ions occupy the A-site, Wyckoff position 2a, with lo- 
cal site symmetry 32, while the magnetic Fe-ion occupies 
a site with inversion symmetry, the 5-site, position 2b, 
site symmetry 3. The Fe-spins order ferromagnetically 
within antiferromagnetically coupled (111) planes with 
the magnetic easy axis perpendicular to [111]. Although 
PE BiFeOs satisfies the magnetic criterion, X transforms 
each Fe sublattice onto itself, XL = T{S\-S2) = L, as can 
be seen in Fig.llfa). In this case of B-site magnetism the 
structural-chemical criterion is not met and the invariant 
^PLM is forbidden by symmetry (the PE point group is 
2'/m' or 2/m for which wFM is allowed). In contrast, 
consider the case where the magnetic ion is on the A-site 
and similarly ordered so that the magnetic criterion is 
still satisfied, Fig.fljb). Now, in this case of A-site mag- 
netism, Xh = -L. Placing a FE- active ion such as Ti^+ 
on the B-site would then satisfy the structural-chemical 
criterion. The PE magnetic point group is now 2/m' 
(2ym) in which wFM is forbidden and by design a FE 
distortion, via £^plm, would lower the symmetry to m' 
(m) thereby inducing wFM. The question remaining is 
whether compounds in our rationalized structure can be 
synthesized? As we discuss next, several already exist. 

The mineral Ilmenite FeTiOa is one member of a 
family of compounds [16] that include the titanantes 
A^+Ti^+Os with A = Mn-Ni. They are ah AFM in- 
sulators with ordering temperatures Tn~40K-100K. At 
atmospheric pressure they form in the Ilmenite structure, 
space group R3. Ilmenite can be thought of as an ordered 
corundum structure. At high-pressure both MnTiOs and 
FeTiOs have been found to form a quenchable metastable 
LiNbOs LBO-phase, space group R3c 13, [3 • Note this 
LBO-phase is structurally isomorphic to BiFeOs except 
the magnetic and FE atom positions are reversed, for ex- 
ample: MnTiOs^BiFeOs implies Mn^Bi and Ti^Fe. 



This is precisely the structural-chemical criterion out- 
lined above. The remaining question is to identify the 
magnetic ground state in the FE phase to determine if 
the magnetic criterion is satisfied. In the remainder of 
this Letter we present a first-principles study of the FE 
and magnetic properties of the LBO-phase of MnTiOs, 
FeTiOs, and NiTiOs- We demonstrate that these are re- 
alizations of the design criteria and provide a novel simple 
picture of how the interplay of D and P leads to electric- 
field control of wFM. 

Method.— We performed density- functional calcula- 
tions using PAW potentials within LSDA+U [19] as im- 
plemented in VASP [13, 1^. The wavefunctions were 
expanded in plane waves up to a cutoff of 500 eV. Inte- 
grals over the Brillouin zone were approximated by sums 
on a 6 X 6 X 6 F-centered /c-point mesh. Phonons were 
calculated using the direct method. Where noted, non- 
collinear calculations with L-S coupling were performed. 
To find appropriate values of on-site Coulomb U and ex- 
change Jh parameters we performed a series of calcula- 
tions to estimate the Curie- Weiss temperature Ocw as a 
function of U for MnTiOs, FeTiOs, and NiTiOs in the 
ground state Ilmenite structure, space group R3 [22|] . For 
all compounds a value of U = 4.5 eV and Jh=0.9 eV was 
found to give a reasonable account of the measured val- 
ues. It should be noted that the presented results do not 
qualitatively change for reasonable variations of U. 

Ferroelectricity.— In the soft-mode theory of ferroelec- 
tricity, the PE to FE transition is associated with the 
softening of a single unstable infrared-active phonon. In 
the R3c^R3c transition, this is a PE phonon polar- 
ized along [111] of symmetry type A2n- We calculated 
the frequencies of these A2n phonons at T=0 and found 
one highly unstable mode, e.g., in MnTiOs, c^so/^ ^1150 
cm~^. The character of this soft mode consists of A-ion 
and Ti displacements moving against oxygen, similar to 



TABLE I: Theoretical (experimental) properties of multifer- 
roic ATiOs, space group R3c. f.u. = formula unit. 





Mn (Exp. [17]) 


Fe (Exp. [18]) 


Ni 




5.127A (5.205A) 


5.05A (5.I2A) 


4.93A 


Ch 


13.63A (I3.70A) 


I3.52A (13.76A) 


13.65A 


Or 


56.4° (56.8°) 


56.2° (56.0°) 


54.7° 


Tfe 


>1000K 


>1000K 


>1000K 


P 


83 /iC/cm^ 


94 /iC/cm^ 


110 /iC/cm^ 


Sew 


-290 K 


-305 K 


-170 K 


Tn 


135 K 


260 K 


100 K 


L 


4.5 liiB/f-u. 


3.6 ijlb/ fM. 


1.6 ijlb/ f'U. 


M 


-0.002 fiB/f.u. 


-0.03 ^is/f-u. 


0.25 iiB/f-u. 


Ji 


-0.5meV 


-l.OmeV 


-1.5meV 


K 


-0.02meV//.iz. 


-0.74 meV//.i^. 


-0.03 meY/ f.u. 


D 


-0.0004 meV//.i^. 


-0.02 meV//.i^. 


-0.35 meV//.i^. 
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other R3c FEs such as BiFeOs and LiNbOs. 

Next we broke inversion symmetry and performed a 
fuh structural optimization within the R3c space group. 
In Table [J the calculated lattice constants at T = 
are shown to be in excellent agreement with those ob- 
served at room temperature for MnTiOs and FeTiOs 
(R3c NiTiOs has not yet been synthesized.) The dis- 
tortion leading from the PE to the FE structure can 
be decomposed entirely in terms of the soft-mode. The 
atomic displacements are almost equal in magnitude to 
those in LiNbOs, meeting Abraham's structural criteria 
for switchable ferroelectricity [23]. Based on Abraham's 
empirically derived formula relating the magnitude of 
these distortions to the FE transition temperature [1^, 
we estimate FE Tc--1500K-2000K. Finally, using the 
modern theory of polarization 0| we calculated a large 
P ^ 80-100 /iC/cm^ comparable to that of BiFeOs- 

Magnetic structure.— Weak ferromagnetism arises as a 
small perturbation - from the relativistic spin-orbit inter- 
action - to a predominately collinear magnetic state, i.e., 
|J|>>|I)|, where J and D are Heisenberg and DM ex- 
change respectively. This difference in energy scales nat- 
urally separates the problem. As such we first identify the 
collinear state that minimizes the spin-interaction energy 
without L-S coupling, i.e., £^h = — Jij^i ■ Sj, by ex- 
tracting the first four nearest neighbor (n.n.) exchange in- 
tegrals, Jn, from total energy calculations [22[. We found 
the state that minimizes En consists of spins aligned 
ferromagnetically within antiferromagnetically coupled 
(111) planes, consistent with the magnetic criterion out- 
lined above. This magnetic state arises due to a strong 
AFM Ji coupling between n.n. spins in adjacent (111) 
planes. The Neel temperature calculated within mean- 
field theory ^ was found to be Tn- lOOK for MnTiOs 
and NiTiOs and Tn- 250K for FeTiOs- 

For a uniaxial crystal the orientation of the global spin 
axis relative to the crystallographic axis is given to low- 
est order by £^ani = ^isin^(^), where is the angle 
between [111] and L. Depending on the sign of the 
spins lie in a plane perpendicular or parallel to [111], 
in which case wFM is allowed or forbidden, respectively. 
To calculate the single-ion anisotropy constant K we first 
performed a self-consistent density-functional calculation 



with collinear spins, without L-S coupling. Then, using 
the charge density and wavefunctions, we performed a se- 
ries of non-selfconsistent calculations with spin-orbit in- 
teraction included for different orientations of the global 
spin axis. We found that the magnetic easy axis lies in 
the plane perpendicular to [111] with K ranging from 
-0.03 meV for MnTiOs and NiTiOs, to -0.7 meV for 
FeTiOs (the much large anisotropy for FeTiOs is asso- 
ciated with the orbital degeneracy). Symmetry then al- 
lows an additional contribution to the total energy which 
can cause the spins to cant: £^dm = Xl^j Dij • (S^ x Sj), 
where Dij is the Dzyaloshinskii vector. Phenomenolog- 
ically this can be described by £^dlm=D-(LxM) where 
L=6'i-6'2 and M=S'i+S'2. Symmetry requires D to point 
along [111], i.e., parallel/antiparallel to P. Since K<0 
requires L to be in the plane perpendicular to [111] and 
subsequently to D, the induced M is perpendicular to 
P. Once the direction of L is fixed, the sign of M that 
minimizes £^dlm is determined by the sign of D. 

Symmetry allows the spins to cant in the FE phase, 
but do they and by how much? To address this question 
we calculated the self consistent spin-density in the pres- 
ence of the spin-orbit interaction for the FE structure 
with for example the polarization pointing down. The 
spins were initialized in a collinear configuration, e.g. 
= (2g/iBS, 0, 0), = (0,0,0), and then allowed to relax 
without any symmetry constraints imposed. For MnTiOs 
the induced moment was rather weak, M = (0, -0.002, 
0) /iB /formula unit (f.u.), although comparable to the 
canonical weak-ferromagnet Fe20s and still measurable. 
The smallness of this result is not too surprising consid- 
ering that the spin-orbit parameter A vanishes for Mn^+ 
in the atomic limit. In contrast, A is relatively large 
for Fe^+ and Ni^+, and correspondingly the induced mo- 
ments increase dramatically; M = (0, -0.03,0) fis/f-U- 
and M =(O.25,0) iib/I-u. for FeTiOs and NiTiOs re- 
spectively j26| . Finally, we can approximate the strength 
of D from the calculated canting angle and Ji. The re- 
sults are summarized in Table HI 

Next we proceed to elucidate the interaction between 
M and P. Similar calculations as we just discussed were 
performed for the PE structure and for the FE structure 
with P in the opposite (symmetry equivalent) direction. 
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FIG. 3: Superexchange pathways, S1-OA-S2 and 5'i-Ob- 
5*2, between nearest neighbor spins. The net Dzyaloshinski- 
Moriya vector, D12, contains two contributions. Da and D^, 
that cancel in the paraelectric phase. 



complex materials. These electrically-controlled switch- 
able magnets provide fertile ground for additional studies 
of how spin and lattice degrees of freedom interact and 
also hold promise for application in magnetic devices. 

Useful discussions with Matthias Bode, Venkat 
Gopalan, Darrell Schlom, and S.K. Streiffer are acknowl- 
edged. Work at the Center for Nanoscale Materials was 
supported by US DOE, Office of Science, Basic Energy 
Sciences under Contract No. DE-AC02-06CH11357. 



again relaxing the spin-density without symmetry con- 
straints. In the former M vanished confirming that the 
FE distortion is required for the observed wFM while in 
the latter M switched sign. These results are consistent 
with our earlier symmetry arguments for a FE inducing 
wFM, i.e., for the invariant Eplm~P-(LxM) [27^ and 
suggests that the direction of M can be switched be- 
tween stable 180° directions by an external electric field 
that would switch P thereby switching M. One way to 
understand this result is to argue that the sign of the 
DM vector D depends on the direction of P, i.e., P cx D. 
At first this may seem puzzling considering the fact that 
Dij is an axial vector. As Ederer and Spaldin pointed 
out [13], P would have to change the sense of oxygen 
rotation in order for D to change sign. In Fig. [2] we com- 
pare the S-O-S bonds in the PE and up/down FE states, 
centering on spin ^2 and its nearest neighbors in adjacent 
(111) planes. A change in chirality of the S-O-S bonds 
is clearly visible due to a change in the direction of P. 
The physics becomes clear by examining the microscopic 
Edm and realizing that the net Dzyaloshinskii vector 
Dij = Dfj has to be summed over all distinct S'l-O^- 
5*2 pathways. In the PE structure two pathways, a left 
and right chiral, contribute to D12 as we show in Fig. [31 
The orientation of the vector is given by 28|, |29[ 
D^2 ^ ^lA X where riA is the unit vector pointing 
along the Si - Oa direction. There is, however, an addi- 
tional pathway connecting to S2 through O^- In the 
PE phase D"^ and D^ can be shown to have equal mag- 
nitude but opposite sign leading to a vanishing net DM 
interaction. In the FE phase, P strengthens/weakens one 
pathway over the other leading to a finite DM interaction 
(in Fig. [2] we only show the "strong" S'i-0-5'2 pathway in 
the FE phase). Therefore the origin of ferroelectrically- 
induced wFM in this class of materials is a change in the 
relative contribution of two DM superexchange pathways 
(with opposite sign) due to a polar lattice distortion. 

Today, the challenge in multiferroics has shifted from 
finding new magnetic ferroelectrics to identifying mate- 
rials in which the polarization and the magnetization are 
strongly coupled. In this work we have presented criteria 
that have the potential to advance the discovery of such 
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